Constraints on Light Dark Matter From Core-Collapse Supernovae 
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We show that hght (~ 1 - 30 MeV) dark matter particles can play a significant role in core-collapse super- 
novae, if they have relatively large annihilation and scattering cross sections, as compared to neutrinos. We find 
that if such particles are lighter than ~ 10 MeV and reproduce the observed dark matter relic density, super- 
novae would cool on a much longer time scale and would emit neutrinos with significantly smaller energies than 
in the standard scenario, in disagreement with observations. This constraint may be avoided, however, in certain 
situations for which the neutrino-dark matter scattering cross sections remain comparatively small. 
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The identity of our Universe's dark matter is one of the most 
interesting questions in modern cosmology. Although a wide 
range of viable particle candidates have been proposed, none 
has been confirmed experimentally. The dark matter candi- 
dates most often studied are weakly-interacting particles with 
masses in the ~ 100 GeV to TeV scale; neutralinos in super- 
symmetric theories being one prominent example. 

Such dark matter particles should not be too light, otherwise 
they could not annihilate sufficiently. Still it is possible to con- 
sider light dark matter (LDM) particles, with the right relic 
abundance to constitute the non-baryonic dark matter of the 
Universe, provided one also introduces new efficient mecha- 
nisms responsible for their annihilations. Such annihilations 
into, most notably, e^e~, could correspond to the exchanges 
of new heavy (e.g. mirror) fermions (in the case of light spin-0 
dark matter particles), or of a new neutral gauge boson U (l||, 
light but very weakly coupled f?], and still leading to rela- 
tively "large" annihilation cross sections. 

The subsequent observation by the INTEGRAL/SPI exper- 
iment of a bright 511 keV 7— ray line from the galactic bulge 
1^] could then be viewed as a sign of the annihilations of 
positrons originating from such hght dark matter particle an- 
nihilations |4, "s"]. These LDM particles, explaining both the 
non-baryonic dark matter of the Universe and the 511 keV 
line, may have spin-^ as well as spin-0 0]. They could even 
potentially improve the agreement between the predicted and 
observed abundances of primordial and ''He, as long as 
LDM particles are not too strongly coupled to neutrinos |7|] . 

Although astrophysical sources such as hypernovae have 
also been proposed as the source of the required positrons [SJ, 
such an origin seems in contradiction with the large extent 
of the 511 keV emission zone |9] (even if the uncertainties on 
the low-energy positron propagation and hypernova rate in the 
bulge are such that it is premature to conclude {13). We shall 
therefore focus on the light dark matter interpretation of the 
511 keV line. (Other exotic particle physics scenarios which 
could generate this emission have been proposed in ifTlll .') 

Given the large rate of positrons produced, smaller dark 
matter masses tend to be preferred, to avoid an excessive pro- 



duction of unobserved 7— rays f^. More specifically, if such 
particles are heavier than 20-30 MeV, internal bremsstrahlung 
(and bremsstrahlung) photons are likely to exceed the ob- 
served number of 7— rays from the galactic bulge 1 12]. And, 
if they were heavier than even 3 MeV, the 7— rays generated 
through the resulting e+e" annihilations might also be incon- 
sistent with observations fl^. 

If MeV-scale dark matter particles do exist, they will be 
thermally generated in the core of collapsing stars. The pres- 
ence of these particles can affect thermal freeze out of weakly- 
interacting neutrinos, depending on their mass, and annihila- 
tion and elastic scattering cross sections. 

Ordinary neutrinos stay in thermal equilibrium through 
weak interactions down to temperatures ~ 2 or 3 MeV during 
the expansion of the Universe, and down to ~ 8 MeV or so, 
in supernovae explosions. Light dark matter particles of mass 
mx annihilate into ordinary ones, staying in equilibrium un- 
til they decouple. This occurs, during the expansion of the 
Universe, at Tp — mx/xp, with xp — 17. In a supernova 
explosion LDM particles will remain in chemical equilibrium 
with other particles, until the temperature drops down to some 
value Tdms. to be determined later (cf. Fig.Q. 

As long as their abundance remains sufficient, these light 
dark matter particles can also influence the behavior of neu- 
trinos in a supernova by having relatively "large" interactions 
with them, e.g. through U exchanges |1] |22]. Neutrinos 
may then be kept longer in thermal equilibrium as a result 
of stronger-than-weak interactions with LDM particles, so 
that their decoupling temperature, in supernovae explosions, 
would be significantly lower than in the Standard Model, if 
dark matter particles are sufficiently light. A crucial ingre- 
dient for this discussion will then be the magnitude of the 
neutrino-LDM elastic scattering cross section. 

We now consider quantitatively these effects through a sim- 
ple model based on the diffusion approximation, largely fol- 
lowing Ref. 1 15]. We begin with the transport equation: 
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where (p is the LDM flux, n the number density of LDM par- 
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tides and ricq its equilibrium value for a LDM mass mx at 
temperature T with vanishing chemical potential. cTann is 
the LDM, anti-LDM annihilation cross section (or self-annihi- 
lation cross section if the LDM is its own antiparticle [23]), 
and Vr is representative of the relative velocity of the two an- 
nihilating particles. 

In the following, cross sections will be taken at typical ther- 
mal energies. We now adopt the diffusion approximation, 
(j) = —D\/n, where D — Aw/3 is the diffusion coefficient 
and A = {J^i ni<Jxi)'^ the LDM mean free path, which de- 
pends on the densities rii of all the particle species (nucleons, 
electrons and positrons, z^'s and i>'s, and in principle dark mat- 
ter particles as well) with which the LDM interacts, with cross 
sections axi- Assuming spherical symmetry and stationarity 
(n = 0), we can write the transport equation as: 
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,Vr (n^ - nl^) , (2) 



where the primes denote derivatives with respect to radius r. 

We define the "LDM-sphere" as the surface beyond which 
LDM annihilations (into e+e~, ...) "freeze out" such that 
the LDM number is effectively conserved further out. The 
number density of LDM particles inside of this LDM-sphere 
should approach its equilibrium value ricq. The radius of the 
LDM-sphere, i?DMS, can be estimated by solving 
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The radius of the surface of last scattering of LDM particles 
is found by solving J^^^dr X)i '''^j — 1- We shall 
concentrate here on LDM scatterings on nucleons (in practice 
mostly neutrons) with density njv(r), and elastic cross section 
axN, the actual Rls radius being at least as large as the one 
we shall estimate by disregarding the other species. This con- 
servative assumption is sufficient to demonstrate that LDMs 
(and therefore eventually neutrinos, with which these LDM 
particles are normally coupled) decouple at lower densities 
and temperatures than in the Standard Model \ 24\. 

For the situations of interest to us, the LDM sphere lies 
within the last scattering sphere, so that the diffusion approx- 
imation is justified. Indeed, the annihilation cross sections of 
LDM particles are normally comparable to LDM scattering 
cross sections with ordinary particles li25ll . so that if LDMs 
can still annihilate they can still also scatter 

To determine the LDM-sphere and surface of last scatter- 
ing for a light dark matter particle, we must adopt a LDM 
mass and a set of (annihilation and scattering) cross sections, 
as well as a distribution of nucleons n n (r) and their tempera- 
ture T(r) in the proto-neutron star. For the latter, we will use 
the following parameterizations, which should be reasonable 
within the range ^15-100 km we are concerned with 1.161 : 



njv(r) ~ < 



6 X 10'''' cm 
6 X lO^'^ cm 



7.8 

' 23km \ , r < 23km, 

12.8 



3 I 23 km 



, 23 km < r < 44 km, 



L2 X 10^3 cm-3 ' 



r > 44 km , 



and 
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5.2MeV(2^) , r<25km, 

1.2 (5) 
5.2MeV(2^) , r>25km, 



r(r) 



These conditions are characteristic for the first few seconds 
over which most of the cooling takes place, in the standard 
scenario. 

We note that the presence of LDM could considerably mod- 
ify profiles compared to the standard scenario. However, we 
will find (cf. Fig.QJ that LDM reproducing the relic den- 
sity and having (Txn ^ o-ann are so strongly coupled that 
they essentially stay in equilibrium as long as they are not 
Boltzmann suppressed, so that they freeze out at temperatures 



DMS 



< 



mx/3 for LDM masses of interest here. 



Annihilation cross sections. The magnitude of the annihi- 
lation cross section of LDM particles into ordinary ones, at 
cosmological freeze out time, is fixed by the relic density re- 
quirement 1 1, 6]. This leads to (cTanni^r/c)? — a few pb. 

In the preferred case of a P-wave annihilation cross sec- 
tion (or at least P-wave dominated at freeze-out), a^nnVr 
is roughly proportional to v^, so that (cranni'r)p-wavc — 
(3 X 10"'^^ cm^) v^/c^. This also turns out to be the right 
order of magnitude for a correct injection rate of positrons 
in the galactic bulge. More precise statements about the re- 
spective roles of P-wave and S'-wave contributions to LDM 
annihilations in the galactic bulge depend, of course, on more 
specific assumptions for the profile of the dark matter mass 
and velocity distribution adopted within the bulge 1 17]. 

We shall often have in mind a simple situation in which a 
scalar or fermionic dark matter particle annihilates (or inter- 
acts) through the virtual production (or exchange) IjJ, |^ of 
a new light gauge boson, U |2], although the present anal- 
ysis is more general. This leads naturally to a P-wave an- 
nihilation cross section, both in the spin-0 case, and in the 
spin-i case as well if the U boson has vectorial (or mostly 
vectorial) couplings to leptons and quarks. Axial couplings 
are already strongly constrained from the non-observation 
of axionlike particles, and parity-violation effects in atomic 
physics 1 2, 6, 18]. For spin-0 particles there may also be S'- 
wave contributions to the annihilation amplitudes, from the 
exchanges of new heavy (e.g. mirror) fermions. 

Results on the temperature of the LDM sphere. In Fig. ^ 
we plot the temperature of the LDM-sphere in the case of 
a P-wave dominated annihilation cross section (oc v^), nor- 
malized to generate the measured dark matter relic density. 
The results for a S'-wave dominated cross section are found 
to be very similar, since it has the same value as a P-wave 
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FIG. 1: Temperature of the LDM-sphere Tdms for cTann^'r oc 
(P-wave), normalized so as to lead to the observed relic density of 
dark matter. The dashed lines represent various hypothesis for the 
LDM-nucleon elastic scattering cross section, taken to be oc 
and (from top to bottom) 1, 10, 10^, 10^, lO", and 10^ times larger 
than the corresponding neutrino-nucleon cross section. The solid line 
shows the case of a massive Vt , for comparison. 

dominated one (up to a factor ~ 2) for a dark matter velocity 
equal to its freeze-out value, v-p ~ 0.4 c. The dashed lines in 
Fig- n correspond to various elastic scattering cross sections 
(see caption). The solid line shows, for comparison, the case 
of a weakly-interacting Vr with a MeV-scale mass. While the 
temperature T ~ 10 MeV resulting for massless neutrinos in 
Fig- n comes out a bit higher than the value ~ 8 MeV from 
more detailed treatments 1 19], what is most important is the 
relative value of the LDM and neutrino temperatures. 

This shows that MeV scale LDMs will remain in equilib- 
rium throughout the proto-neutron star at least down to rel- 
atively low temperatures T ~ 3 MeV, as an effect of the 
large values of the annihilation cross sections of LDM par- 
ticles (into ordinary ones). This occurs even if we do not 
assume rather high values of the scattering cross sections of 
LDM particles with ordinary ones. Large scattering cross sec- 
tions then contribute to further reinforce the effect by increas- 
ing the LDM diffusion time allowing to keep LDM particles 
at chemical equilibrium down to even lower values of the tem- 
perature, possibly down to Tdms — 1 MeV, as illustrated by 
the lower dashed curves of Fig.Q] 

Consequences for the neutrino temperature. Thus light 
dark matter particles with relatively large annihilation cross 
sections (as required from relic abundance) remain in equi- 
librium down to lower temperatures, T < 3 MeV. This fea- 
ture may be transmitted to neutrinos, that will themselves stay 
longer in thermal equilibrium as a result of their interactions 
with LDM particles, provided neutrino-LDM cross sections 
are also enhanced as compared to ordinary neutrino cross 
sections. 

The kinetic equations ( I2I3> and the one fixing R^s are for- 



mally the same for neutrinos, substituting the relevant cross 
sections for scattering and annihilation of neutrinos and the 
equilibrium density of the relevant neutrino flavor. Inside the 
LDM sphere the relevant quantities for neutrinos may be ap- 
proximated as 



3 [nNCJvN + ?^oq,x cFvx) 

lyiJ^XX [or XX) + '''SM 



(6) 



where cross sections are for the processes indicated as sub- 
script (taking also into account the anti-LDM contribution if 
LDM particles are not self-conjugate), and ctsm indicates the 
Standard Model contribution. Note that the LDMs are kine- 
matically accessible by neutrinos for temperatures not much 
lower than mx /3. If indeed the cross sections for neutrino- 
LDM scattering and neutrino annihilation into LDMs are 
comparable to the ones for LDM-nucleon scattering and LDM 
annihilations into leptons (supposed to be "large"), respec- 
tively, the quantities in Eq. (|6j will be dominated by the non- 
standard contributions. This is because LDM cross sections 
are normally ~ a few pb (at freeze-out velocity, to give the 
appropriate relic abundance), more than ~ 10* larger than 
weak-interaction cross sections ( ~ G^p ), at the relevant 
energies. 

This implies that neutrinos (if indeed they have relatively 
"large" interactions with LDM particles) should stay in chem- 
ical equilibrium at least as long as the LDMs do and T > 
mx/'i. We then conclude that mx ^ 10 MeV would give 
rise to neutrino decoupling temperatures < 3.3 MeV for all 
flavors, as compared to ~ 8 MeV for and Vr in the stan- 
dard scenario. 

This would make it quite unlikely to observe neutrinos with 
energy of order 30 — 40 MeV, as have been observed from 
SN1987A 1 20], especially for emission spectra that are sup- 
pressed at the highest energies compared to thermal distribu- 
tions because the cross sections increase with energy 1 19], in 
which case we can conclude that lighter LDM masses < 10 
MeV are practically excluded. 

All this relies, of course, on the potentially "large" size 
of the neutrino-LDM scattering and vv LDM's annihila- 
tion cross sections, normally expected to be comparable to 
the "large" LDM's e+e^ annihilation cross section. 

It is worth noting, however, that there are special situations 
for which the U boson would have no coupling at all (or sup- 
pressed couplings) to neutrinos i2lll . And that for a spin-0 
LDM particle interacting with ordinary ones through the ex- 
changes of heavy (e.g. mirror) fermions, the i^-LDM interac- 
tions would be severely suppressed (as compared to electron 
or nucleon-LDM interactions), as a result of the chiral charac- 
ter of the neutrino field 1 1, 6]. In both cases we end up with 
no significant enhancement of neutrino-LDM interactions, so 
that the presence of the LDM particles has no direct signifi- 
cant effect on the behavior of neutrinos, then still expected to 
decouple at ~ 8 MeV (for i/^ and ;/,-), as usual. In such a 
case, no new constraint is obtained on the mass mx of LDM 
particles. 
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Furthermore, the above resuhs may also be obtained, or un- 
derstood, as follows. Let us return to LDM particles rather 
"strongly" coupled to neutrinos (and nucleons), both types of 
particles decoupling at T < 3.3 MeV. As LDMs can then 
contribute, at most, as much to the cooling flux as the neu- 
trinos (due to fewer degrees of freedom), the cooling time 
scale would be larger than in the standard scenario by a factor 
> (8/3.3)'^/2 ~ 20 because the thermal flux is also oc T'^. 
As SN1987A observations were consistent with the standard 
cooling time scale of 10-20 s, such non-standard scenarios are 
then very strongly disfavored, to say the least. 

The cooling time scale can also be estimated by the dif- 
fusion time Tdiff ^ R^g/X. This is dominated by the 
innermost regions of the hot neutron star of size i?NS — 
10 km, whose density is not significantly modified by the 
presence of LDM. At a typical temperature T ~ 30 MeV, 
a^^N - llG|rV7r ^ 1-7 x 10""^" cm^, and at nucleai- 
densities Ucq ^ npf/100. Thus, for neutrino-LDM cross 
sections comparable to electron-LDM ones (i.e. typically > 
4 X lO^'^^cm^, so that a,,x ^ a,^N), the neutrino mean 
free path is dominated by interactions with the thermal pop- 
ulation of LDMs, so that Xi, ^ (ricqcri^x)^^ ^ 0.3 cm, 
as compared to X„ ~ inNa,^N)~^ ^ 35 cm in the stan- 
dard scenario. The LDM mean free path is even shorter, 
Aldm ^ {nNcrxN)~^ ^ 1-5 X lO"'^ cm, assuming ctxn ^ 
Ci/X ^ cTeX (or even less if LDM self-interactions were to 
contribute significantly). In the interior of the proto-neutron 
star the energy flux is thus dominated by neutrinos. The cool- 
ing time scale is a factor > 100 larger than in the standard 
scenario, consistent with the previous argument. This cooling 
time argument may be extended up to higher LDM masses 
~ 20 or even 30 MeV, i.e. as long as LDMs are significantly 
present at T ~ 30 MeV, and rather "strongly" coupled to neu- 
trinos. 

Given that about 3 x 10^'^ erg of binding energy has to 
be liberated during Tdiff , in the relativistic regime the freeze 
out temperature will scale as Ti, oc t^^'^. For gxn ^ 
lO"' a^x {T ~ 30 MeV), this argument suggests will be a 
factor > 3 times smaller than usual, as found previously. 

In conclusion, we have demonstrated that light dark mat- 
ter models with generically "large" cross sections fixed by re- 
quiring them to reproduce the relic dark matter density are 
considerably disfavored by the resulting modification of core- 
collapse supernova cooling dynamics if the dark matter mass 
is < 10 MeV, at least. 

Depending on how strict 7— ray constraints from the galac- 
tic bulge are, the new supernovae constraint presented here 
could strongly disfavor the possibility that annihilating dark 
matter particles be the source of the 511 keV emission from 
the galactic bulge. 

Or, conversely, these new results could indicate that 
neutrino-LDM interactions should not be enhanced, favoring 
a U boson with no (or small) couplings to neutrinos and/or a 
spin-0 dark matter particle interacting through heavy fermion 
exchanges. 
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